Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5th February 2008 (MN WT^ stylo file v2.2) 



Kinematics of the ultra- compact helium accretor AM 
Canum Venaticorum 

G.H. A. Roelofs/* P.J. Groot,^ G. Nelemans,^ T.R. Marsh^ and D. Steeghs' 

- - - ^Department of Astrophysics, Radboud University, Toernooiveld 1, 6525 ED Nijmegen, The Netherlands 
Department of Physics, University of Warwick, Coventry CV4 7AL, UK 
^ Harvard- Smithsonian Center for Astrophysics, 60 Carden Street, Cambridge, MA 02318, USA 

o 

(N 

■ 5th February 2008 

=3 

1—5 

Cn ' ABSTRACT 

We report on the results from a five-night campaign of high-speed spectroscopy of the 
17-minute binary AM Canum Venaticorum, obtained with the 4.2-m Wilham Herschel 
Telescope on La Palma. 

We detect a kinematic feature that appears to be entirely analogous to the 'central 
spike' known from the long-period, emission-line AM CVn stars GP Com, V396 Hya 
and SDSS J124058. 03-015919. 2, which has been attributed to the accreting white 
. dwarf. Assuming that the feature indeed represents the projected velocity amplitude 

' and phase of the accreting white dwarf we derive a mass ratio q = 0.18 ± 0.01 for 

I AM CVn. This is significantly higher than the value found in previous, less direct 

■ measurements. We discuss the implications for AM CVn's evolutionary history and 

show that a helium star progenitor scenario is strongly favoured. We further discuss 
the implications for the interpretation of AM CVn's superhump behaviour, and for 
the detectability of its gravitational-wave signal with LISA. 

O I In addition we demonstrate a method for measuring the circularity or eccentricity 

of AM CVn's accretion disc, using stroboscopic Doppler tomography. We test the 
^ ■ predictions of an eccentric, precessing disc that are based on AM CVn's observed 

' superhump behaviour. We limit the effective eccentricity in the outermost part of 

the disc, where the resonances that drive the eccentricity are thought to occur, to 
' e — 0.04 ± 0.01, which is smaller than previous models indicated. 

H I Key virords: stars: individual: AM CVn - binaries: close - novae, cataclysmic vari- 

ables - accretion, accretion discs 
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1 INTRODUCTION 

AM Canu m Venaticorum was fo u nd in a survey of faint, blue 
objects by iHumason fc Zwickvl l^i^)- It was observed to 
have peculiarly broad and shallow helium absorption lines, 
but no hydrogen (Grccnstcin 1957). The star was sh own to 
be a possible ultra-compact white dwarf binary by ISmakI 
il967l) . who discovered photomet ric variations on an 18- 
minute period; quickly thereafter, IPaczynsTil lll967f) noted 
that it would be a prime example of a binary whose evolution 
is expected to be governed by gravitational-wave radiation, 
and which could serve as an excellent test for the existence 
of gravitational radiation. The interpretation as an interact- 
ing binary a nalogous to the Cataclysmic Variables was first 
proposed by 'War ner fc RobinsonI I^T^) upon their discov- 
ery of rapid flickering in the light-curve. The discussion as 
to the true orbital period of the system was finally put to 
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rest by iNelemans. Steeehs fc Grood i200ll) who discovered 
a kinematic 'S-wave' feature in time-resolved spectra of AM 
CVn, thereby proving that the orbital period is 1028 sec- 
onds, while the main photometric signal at 1051 seconds is 
to be interpreted as a 'superhump'. 



Here we present phase-resolved spectra of AM 
CVn with significantly higher spectral resolution and 
signal-to-noise ratio than thos e previously used by 
INelemans. Steeehs fc Grood (1200 J) . Moreover, our current 
data-set fully samples the proposed 13.37-h accretion disc 
precession period (Patterson. Halpern fc Sh ambrook 1993; 
ISkillman et'al]ll99if for the first time, allowing for a char- 
acterisation of the spectroscopic appearance of the system 
as a function of orbital period, precession period, and super- 
hump (beat) period. 



2 G.H.A. Roelofs et al. 



Date 


UT 




Exposures 
(30-40 s) 


Typical 
seeing (") 


2005/04/20 


20:52- 


-05:06 


523 


0.8 


2005/04/21 


22:46- 


-05:15 


114 


0.8 


2005/04/22 


20:41- 


-05:04 


585 


1.1 


2005/04/23 


21:03- 


-04:57 


370 


1.4 


2005/04/24 


20:35- 


-05:01 


567 


0.8 



Table 1. Summary of our observations. 



Line Equivalent width (A) 

Hei 4387 1.58 ± 0.02 

He I 4471 {-l-He I 4437) 2.84 ± 0.02 

Hen 4686 (+Hoi 4713) 0.76 ±0.04 

He I 4921 2.60 ±0.05 

He I 5015 0.73 ±0.02 



Table 2. Equivalent widths of the detected accretion disc absorp- 
tion lines, including estimated errors, in the average spectrum of 
AM CVn (figure 



2 OBSERVATIONS AND DATA REDUCTION 

We obtained phase-resolved spectroscopy of AM CVn on 
20, 21, 22, 23 and 24 March 2005 with the William Herschel 
Telescope (WHT) and the ISIS spectrograph. The observa- 
tions consist of 2159 spectra taken with the 1200B grating, 
covering ~4300-5100A; the exposure time was 30-40 sec- 
onds depending on sky transparency and airmass. The sky 
varied from clear to opaque due to high clouds, but the me- 
dian seeing was quite good, about 0.8", giving an effective 
resolution of about 0.7 A or 45 km/s. 

The chosen slit width was dependent on the seeing and 
transparency, and varied from 1.0" — 1.4". The detector was 
the standard EEV12 chip in the spectrograph's blue arm, 
windowed and binned by a factor of two in the spatial di- 
rection to increase the read-out speed. The read-out speed 
was kept 'low' to minimize the read-out noise. Each night, 
an average bias frame was created from 20 individual bias 
exposures, and a normalised flatfield frame was constructed 
from 20 incandescent lamp flatfields. 

All spectra were extracted using the IRAF implementa- 
tion of optimal (variance-weighted) extraction. The read-out 
noise and photon gain, necessary for the extraction, were cal- 
culated from the bias and flatfield frames, respectively. He- 
NeAr arc exposures were taken every hour during the night 
to correct for instrumental flexure; all dispersion solutions 
were interpolated between the two arc exposures nearest in 
time. Arcs were taken before and after every rotator adjust- 
ment needed to keep the slit at near-parallactic angle. In 
each arc exposure a total of about 40 arc lines could be fit- 
ted well with a Legendre polynomial of order 4 and 0.04 A 
root-mean-square residuals. The root-mean-square residuals 
due to flexure are estimated to be smaller than this; the to- 
tal arc drift was measured to be 0.5 A peak-to-peak over an 
entire night. All spectra were transformed to the heliocentric 
rest-frame prior to analysis. 

The average spectrum was corrected for instrumental 
response using spectrophotometric standard star Feige 34. 

A summary of all observations is given in table Q 



3 RESULTS 

3.1 Average and phase-binned spectra 

The grand-average spectrum of AM CVn is shown in figure 
It shows several broad neutral helium absorption lines 
seen before in AM CVn's spectrum, and clear emission of 
Hell 4686. The absorption wings extend to about 1300 km/s 
(half width at half minimum). See table|5|for a list of equiv- 
alent widths for the strongest absorption lines. 




4400 4600 4800 

Wavelength (A) 



Figure 1. Average spectrum of AM CVn. The wavelength- 
dependent instrumental response has been corrected for, but the 
overall received flux is strongly affected by clouds; no attempt has 
been made to correct for this. 

The phase-binned (trailed) spectrum reveals several ad- 
ditional, weaker features that are too smeared out in the av- 
erage spectrum to be detected. See figures |5| and |21 Among 
these are the 'S-wavo' emission features in several lines, al- 
ready found by Nclcmans, Stccghs & Groot (2001). It in- 
cludes the m etal line Mgll 4481; already prese nt but not 
mentioned in lNelemans. Steeehs fc Grood i200ll) . Two very 
weak S-wave features flank the (also very weak) He I 5048 
line, which we identify with Sill 5041 & 5056. 

The most remarkable is a weak and narrow emission 
feature slightly offset to the red of He I 4471. It is clearly 



Hel \ 4387 Hel A 4921 
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Figure 2. Phase- binned, average-subtracted spectrum of AM 
CVn around the Hel 4387 and Hel 4921 lines. 
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ofTset in phase from the presumed bright spot feature, and 
has a much lower amphtude (see figure |^ . As such it has 
everything in common with the 'central spike' feature that 
was first detected in GP Com l,Smak f975.,l and has since 
become a common (though poorly under stood) feature in 
long-period, emissi on-line AM CVn stars (|R,uiz et alJl2nn]l 
iRoelofs et al]|2005^ . As in the aforementioned systems, the 
phase of the low-velocity component (relative to the bright 
spot) is consistent with the expected phase of the accret- 
ing white dwarf, and its relatively low velocity amplitude is 
difficult to reconcile with any emission site in the accretion 
disc. It is also intrinsical ly redshifted through some as yet 
unkn own mechanism (cf. 'Marsh"l99^;| Morales-Rueda et alJ 
l2003l) . Based on these three properties, we interpret the fea- 
ture as the perfect analogue of the central spike in GP Com, 
V396 Hya, and SDSS J124058.03-015919.2. 

We see no evidence for a central spike in any of the ot her 
helium lines. However, from lMorales-Rueda et al] ll2003l) we 
see that the central spike is strongest in the He I 4471 line 
in GP Com, and it is clear that we are barely detecting this 
line's spike in AM CVn. Our non-detection of a spike in the 
other helium lines is thus plausibly a matter of signal-to- 
noise. 



3.2 Doppler tomography 

3.2.1 The central spike 

The weak, sinusoidal emission feature we interpret as the 
central spike in the He I 4471 line can be made more 
tangible in a linear back-projection Doppler tomogram 
jMarsh fc Horn(^1988^ . We begin by phase-binning our 2159 
spectra in ^^ 200 bins using the known orbital period P prb = 
1028.7322 s iHarvev et alJll99d ISkillman et al.lll999^ ■ The 
two bright spots close to the spike - from the He I 4471 
line itself and from the nearby Mgll 4481 line, see figure|21- 
are masked out to prevent them from dominating the (back- 
projected) spectra. Each phase-binned spectrum is then nor- 
malised by fitting a polynomial and dividing the spectrum 
by it, and finally all phase bins are divided by the column- 
averaged spectrum so that we are left with the components 
of the spectrum that vary on the orbital period, apart from 
the bright spots which were masked out earlier. 

The next step is determining the central spike's rest 
wavelength or, equivalently, its intrinsic (rest) velocity 7 rel- 
ative to the He I 4 471 rest wavelength. We proceed in the 
same manner as in iRoelofs et al.l ll2005l) . that is, we make 
Doppler tomograms for a range of trial wavelengths and de- 
termine th e wavelength at whic h the spectra auto-correlate 
best. As in lRoelofs et alJ (1200511 . the height of the presumed 
central spike nicely peaks at a certain wavelength, which we 
take as the rest wavelength of the central spike. See table |3 
and the corresponding Doppler tomogram in figure |1| 

The final step is determining the central spike's velocity 
amplitude and phase, and their as sociated error s . We em- 
ploy the bootstrap method used in lRoelofs et alJ i200El) . In 
this Monte Carlo process we make a large number (10^) of 
Doppler tomograms with the recipe described above, using 
a random selection of 2159 spectra taken from our data set, 
and allowing for replacement. We determine the emission 
centroid of the central spike in each tomogram, make a 2-D 
histogram of all the Kx and Ky values obtained, and fit a 



Feature 


7 


Ky 


Kx 




(km/s) 


(km/s) 


(km/s) 


Bright spot 


0±5 


330 ±8 


-336 ± 8 


Central spike 


40 ±4 


-92 ±4 


0±4 



Table 3. Velocities of the bright spot and central spike compo- 
nents of the He I 4471 line, after aligning the central spike with 
the negative Ky axis. 



2-D Gaussian to this distribution. The width of this Gaus- 
sian is taken as a measure for the error on the central spike's 
coordinates in Kx, Ky space. See table|H]for the results. 

3.2.2 The bright spot 

We follow the same procedure as for the central spike, but 
now masking out only the bright spot of the Mgll 4481 line. 
We again construct an ensemble of ~10^ Doppler tomograms 
via the bootstrap method in order to be able to estimate 
the errors on the resulting data points. Table |3| lists the 
rest wavelengths, velocity amplitudes and phases as well as 
their errors, of both the bright spot and the central spike of 
the He I 4471 line. The central spike has been aligned with 
the negative Ky axis, which coincides with the conventional 
zero-phase of the accreting white dwarf. 

Figure shows the velocities of the central spike and 
the bright spot, with errors, to show their relative phases 
and amplitudes, after aligning the central spike with the 
negative Ky axis. Looking ahead at the next section, we 
also plot the velocities of a free-falling accretion stream, as 
well as the Keplerian velocities along the trajectory of the 
ballistic accretion stream, for a binary of mass ratio q = 
M2/M1 = 0.18. 

3.3 The mass ratio of AM CVn 

Using the velocity amplitudes and phases of the central spike 
and the bright spot, we can constrain the mass ratio and the 
effective accretion disc radius in AM CVn. To this end we 
solve the equation of motion for a free-falling stream of mat- 
te r through the inner La grangian point, based on the results 
of lLubow fc Siiul lll975r) . and we see whether the resulting 
accretion stream and/or accretion disc velocities and phases 
at the stream-disc impact point match with the measured 
values. The results are displayed in figure |n] 

Keplerian velocities in the bright spot imply very small 
effective accretion disc radii {R ~ 0.3-Rli, with Rli the dis- 
tance from the centre of the accretor to the inner Lagrange 
point), and very large mass ratios {q ~ 0.5). In particu- 
lar, the disc would not come anywhere near the 3 : 1 reso- 
nance radius. This resonance is commonly thought to drive 
the superhump behaviour that is observed in outbursting 
dwarf novae as well as nova-like systems such a s AM CVn 
JWhitehurst fc Kinelll99ll iHirose fc Osakilll990t) . 

For purely ballistic stream velocities in the bright spot, 
the best-fitting accretion disc radius is 5 — 10% smaller than 
the maximum accretion disc radius that can be contained 
within the primary star's Roche lobe. In this case, the disc 
would clearly extend past the 3 : 1 resonance radius. If we 
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Line velocity (km s ) 

Figure 3. Phase-binned, average-subtracted spectrum of AM 
CVn around the He I 4471 line. The bright spot trails around 
zero velocity (component a), in phase with the bright spot from 
the nearby Mgll 4481 line (b). In addition there is a weak com- 
ponent moving around 40 km/s with a semi-amplitude of 92±5 
km/s, and offset in phase from the bright spot (c). We propose 
that this is the 'central spike' of the He I 4471 line. 



assume purely ballistic stream velocities, we obtain a mass 
ratio q = 0.18 ± 0.01. This constitutes a lower limit. 

If we allow for mixing between the ballistic stream and 
the Keplerian disc velocities, the mass ratio could be slightly 
larger than this. See figure |S| The maximum mass ratio for 
which the bright spot would remain outside the 3 : 1 reso- 
nance radius is g ^ 0.22±0.01, which corresponds to a (best- 
fitting) mix of about 80% stream and 20% disc velocities. A 
mass ratio oi q < 0.25 is commonly quoted as the require- 
ment for a binary to be able to excite its 3 : 1 resonance and 
develop the corresponding superhum ps ( Whi tehursd llQsd 
IWhitehurst fc KindlTgOlt iHirose fc Osaki,.199cE 

Two values for the mass ratio q from the literature 
are overplotted in figure |S| Both these values are de- 
rived from the fractional superhump period excess e = 
(Psh - Pqa) /{orb = 0.0216 observed in AM CVn (e.g. 
ISkillman etai]ll999D . combined with theoretical-numerical 
or empirical relations between e and q. The mass ratio 
derived from AM CVn's kinematics is significantly larger 
than these predictions. The present data indicate a lower 
hmit oi q = 0.18 ± 0.01, compared with q — 0.10 as ob- 
tained from the latest empiri cal e(q) relation for hy drogen- 
rich Cataclysmic Variables bvlPatterson et ah| l| 2005[) fPOS). 
The value q = 0.087 found b\ TNelemans^teegh^g Grood 




K>, (km s ') 

Figure 4. Linear back-projection Doppler tomogram of the red- 
shifted, low-velocity central spike. Its velocity amplitude comes 
out at 92it5 km/s. Its phase has been aligned with the negative 
Ky axis, the conventional zero-phase of the accreting star. 




o 
o 
Ln 



-500 







Kv (km/s) 



500 



Figure 5. The measured central spike and bright spot velocities, 
with the associated 3-cr error bars. Overplotted is the best-fitting 
binary of mass ratio q = 0.18, with the ballistic accretion stream 
starting at the inner Lagrange point and passing through the 
bright spot. The upper solid line starting within the secondary 
star represents the Keplerian accretion disc velocities along the 
ballistic stream trajectory. The labels represent the accretion disc 
radius as a fraction of Rli- 
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l|200ll) (NSGOl), from fits to numeric al accretion disc sim- 
ulations as taken from IWarneJ (Il995ll . deviates even more. 
In the remainder of this paper, we shall conservatively use 
q — 0.18 ±0.01, which corresponds to pure stream velocities 
in the bright spot. 



3.4 Stroboscopic Doppler tomography 

The permanent superhump phenomenon, exhibited by AM 
CVn, is usually explained in terms of a precession of a tidally 
deformed (i.e., non-circular) accretion disc. After exactly one 
orbital cycle, the orientation of the two stars with respect to 
the observer is restored, but during the orbit the eccentric 
disc has advanced slightly so that it takes slightly more than 
one orbital period to restore the orientation of the two stars 
with respect to the accretion disc. This beat period is then 
assumed to be responsible for the main photometric signal, 
which is usually observed to be a few percent longer than 
the orbital period. 

The model requires that there be a precessing ac- 
cretion disc with s ustained eccentricity. Furthermore, 
ISkillman et alJ 1^9^, among others, find the superhump 
cycle to be extremely stable over many thousands of orbital 
periods, with the superhump period always within 0.1 s of 
1051.2 s and the superhump waveform remaining very sta- 
ble throughout the years^ . This implies that the shape of the 
disc is (semi-)stationary rather than chaotic. It further im- 
plies that the superhump phase will drift by at most 0.04 of 
a cycle over our ~100-hour observing baseline, if we adopt 
their average superhump period of 1051.2 seconds. These 
two results open up the possibility of combining the large 
number of superhump cycles in our data-set to reconstruct 
the appearance of the system as a function of superhump 
phase. 

The application of Doppler tomography as an instru- 
ment for testing the eccentric disc model is as follows. As 
the secondary star goes around the slowly precessing and 
non-circular disc, the effective radius at which the accretion 
stream impacts the disc will vary. Since the bright spot is 
caused by the free-falling stream of matter from the inner 
Lagrange point crashing into the edge of the accretion disc, 
a variation in effective accretion disc radius will lead to dif- 
ferences in the gas velocities of both the stream and the disc 
matter at the impact point. Observationally, we then ex- 
pect the bright spot to show excursions in velocity space, as 
a function of the relative orientation of the eccentric disc and 
the secondary star, i.e., as a function of superhump phase. 

In order to measure these excursions we must divide 
our spectra in a number of bins corresponding to difi'erent 
superhump phases, and subsequently phase-fold the spectra 
contained in each bin on the orbital period. The method 
bears resemblance to t he 'stroboscopic Dopp ler tomography' 
method employed by iMarsh fc Duck! lll996ll to disentangle 
the orbital and white dwarf spin periods in the magnetic 
accretor FO Aqr, except that we now have the orbital period 
of the binary and the precession period of the accretion disc. 

We argue that this is in fact quite a sensitive method 



^ Note that the superhump waveform is double- humped; most of 
the power is actually in the first harmonic at 525.6 seconds. 



for measuring variations in the effective radius of the accre- 
tion disc. If we consider the simplest case of an eccentric, 
elliptical accretion disc, a simple estimate for the ratio of 
the projected apastron and periastron velocity amplitudes 
"f^a.p as a function of the eccentricity e is 

!^^i±£^l + 2e, e«l. (1) 

Va, 1 - e 

A bright spot of average velocity amplitude v ~ 450 km/s, 
combined with a reasonable spectral resolution and signal- 
to-noise ratio that allow the bright spot to be measured to 10 
km/s accuracy, gives an eccentricity sensitivity of e ~ 0.01. 
In practice, the excursion of the bright spot in a Doppler 
tomogram will be larger and the eccentricity sensitivity will 
correspondingly be better than this. 

As an example we show in figure Q by solving numer- 
ically the equation of motion for a ballistic test mass in 
the binary potential, how we expect the bright spot to vary 
as a function of superhump phase, in a binary of mass ra- 
tio q = 0.18, with an elliptical accretion disc of semi-major 
axis R = 0.65-Rli and eccentricity e = 0.2. This eccentric- 
it y value is chosen based on the simp le observational model 
of|PattCTSon^HalDern & Shambrook and the numer- 

ical models of^impson & Wood (1998). The bright spot is 
expected to show excursions of several hundreds of km/s, 
indicating that eccentricities of this magnitude should be 
easily observable with our spectral resolution. 

In our simple model we have assumed a somewhat sharp 
edge to the accretion disc as in e.g.^rmita gc fc Livio ( 1998), 
such that the gradual changes in pressure along the ballistic 
stream do not significantly alter the depth of the stream's 
impact into the disc. In order to explore the possible ef- 
fects of such variations, consider the pressure of the accre- 
tion stream ps, 

1 2 

Ps — Pstatic ± Pdynamic ^ Pdyiiamic — 

where p and v are the density and velocity of a fluid element 
along the stream, respectively, and Pstatic can be neglected 
since fcsT ^ rriYicV^ in the free-fall regime (by definition, 
one could say), where v is of the order 10^ km/s and the 
temperature T of the stream, set by the surface temperature 
of the donor star, is of the order 10* K. If we assume a steady 
ballistic fiow of gas, where the scale height of matter in the 
flow rem ains fairly consta nt in the region just outside the 
disc (.Lubow fc Shulll976D since the matter does not have 
time to maintain hydrostatic equilibrium perpendicular to 
the flow, then conservation of mass gives 

p oc v"^ (along streamline) (3) 

so that 

Ps oc u (along streamline) (4) 

The pressure along the accretion stream thus increases 
linearly with velocity. In case of a 'soft' edge to the accretion 
disc (that is, a slowly increasing disc pressure Pd from the 
disc edge inwards) this could have the effect of increasing 
the ratio of observed periastron and apastron bright spot 
velocities for given e, as the stream impacts slightly deeper 
into the disc rim at periastron due to its higher pressure 
than it does at apastron. Thus if one employs the radius at 
which the bright spot occurs as the effective in situ accre- 
tion disc radius, and one allows the disc edge to be soft, an 
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q = M2/M, q = Mj/M, q = Mj/M, 

Figure 6. Mass ratio q and effective accretion disc radius R of AM CVn based on the observed velocities and phases of the central 
spike and the bright spot. Left: purely ballistic stream velocities; right: purely Keplerian disc velocities; centre: mixing of 80% stream 
and 20% disc velocities in the bright spot. The grey-scale shows the exclusion confidence levels, in standard deviations. The upper and 
lower dotted lines indicate the edge of the primary Roche lobe and the circularisation radius, respectively, while the dashed line shows 
the 3 : 1 resonance radius. The labels NSGOl and P05 in dicate previous mass ratio estimates based on AM CVn's superhump behaviour 
jNelemans. Steeghs fc Groo jEoOll: IPatterson et alJl200Rft . 
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Figure 7. Simulated stroboscopic Doppler tomogram for a bi- 
nary of mass ratio q = 0.18, accretion disc semi- major axis 
R = 0.65i?Lli and accretion disc eccentricity e = 0.2. The ellip- 
tical group of spots represent Keplerian accretion disc velocities 
in the bright spot, the group sticking out to the left represent 
ballistic accretion stream velocities, as a function of superhump 
phase (one full cycle shown). 



eccentricity measurement based on the bright spot excursion 
amplitude could lead to an overestimate of e, and thus to 
an upper limit for e. 



3.5 Eccentricity of the accretion disc of AM CVn: 
the bright spot 

Figures |S| and show the results from our stroboscopic 
Doppler tomograms. In fig. |S|we plot eight Doppler tomo- 
grams that together cover one full superhump period. In 
order to improve on the signal-to-noise ratio of the individ- 
ual maps we combined the data (in velocity space) of four 
lines that show a clear bright spot in the trailed spectra, 
namely He I 4387, He I 4471, Mgll 4481 and He I 4921 (fig- 
ures Figure |5| shows the measured excursions of the 
bright spot as a function of superhump phase, for 16 super- 
hump phase bins, connected by a line, where the bright spot 
centroids have been determined in the same way as in sec- 
tion 13.21 Here, the 16 phase bins oversample the superhump 
cycle by a factor of two, so that each bin contains l/8th of 
the data as in figure |S| 

Based on the observed bright spot excursion amplitude 
of less than 100 km/s (see figure|5J, we constrain the effective 
eccentricity of the outermost part of the accretion disc to 
e = 0.04 ± 0.01, where we define the 'efi'ective' eccentricity 
as that which would cause a simple elliptical disc to give 
rise to the same bright spot excursion amplitude in velocity 
space. The error on this value can be rather small because 
the method of measurement is insensitive to the accretion 
disc's semi-major axis (or average radius), and to whether 
the bright spot represents stream or disc velocities (see also 
figure|7|l. In case the bright spot velocity is a variable mix of 
the disc and stream velocities, the derived eccentricity can 
be an overestimate. 

Despite the rather large error bars in the left panel of 
figurel^l at least relative to the small movement of the bright 
spot, the disc appears to be more complex than the simple 
ellipse which it is often pictured to be in explaining the su- 
pe rhump phenomenon. Deta iled numerical simulations, e.g. 
by ISimpson fc WoodI il998l) , show that the accretion disc 
in a permanent superhump system like AM CVn is in fact 
most likely of irregular shape. Not only that, the disc is also 
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Figure 8. Doppler tomograms showing the bright spot as a func- 
tion of superhump phase. One full cycle is shown from left to right, 
top to bottom. The bright spot is visible at all superhump phases, 
and shows no large excursions in velocity space. 



predicted to change shape during a superhump period, only 
to return to its original shape after one full superhump cycle 
- therefore, the disc could m principle be exhibiting strong 
changes in radius at the other side of the accreting star, 
where the mass stream cannot be used to probe its effective 
radius. The changes in effective radius of the accretion disc 
might thus be larger than is revealed by the excursions of 
the bright spot. 



3.6 Eccentricity of the accretion disc of AM CVn: 
integrated disc profile 

Figure ^1 shows trailed spectra and Doppler tomograms of 
Hell 4686 as a function of superhump phase. The Hell 4686 



line was chosen because it clearly shows emission in the aver- 
age spectrum (fig.0, which is commonly used as a tracer of 
locally enhanced dissipation in accretion discs, for instance 
in searches for spiral shocks in outbursting dwarf novae. A 
sinusoid with the phase and amplitude of the central spike - 
i.e. the presumed motion of the accretor - has been shifted 
out of the spectra to leave the net kinematic signal of the 
matter in the disc orbiting the accretor. The ionised helium 
emission profile is clearly variable on the superhump period: 
the emission is asymmetric and the degree of asymmetry 
is changing (see also fi gure El. These g eneral features agree 
with the predictions of ISimpson fc Woo d (1998). 

It can furthermore be deduced that the He ll 4686 emis- 
sion is strongest where the effective radius of the disc is 
smallest. The right panel in figure |5| shows both the bright 
spot velocity amplitude and the angle between the bright 
spot and the centroid of He ll 4686 emission, as a function of 
superhump phase. This angle is given in spatial coordinates 
- for the transformation of velocity into spatial coordinates 
we assume roughly ballistic stream velocities for the bright 
spot based on the results of section and roughly Keple- 
rian disc velocities for the Hen 4686 emission. Despite the 
rather large error bars, there appears to be a trend whereby 
the maximum bright spot velocity amplitude (which should 
correspond to the smallest effective disc radius) corresponds 
to alignment of the bright spot with the He ll emission cen- 
troid, while the minimum bright spot velocity (i.e. largest 
effective disc radius) corresponds to anti-alignment. 

These findings again ag r ee wi th the numerical simu- 
lations of ISimpson fc Woodl il998h . Their analysis of the 
changes in internal energy in the disc as a function of super- 
hump cycle, too, showed that energy production was highest 
where the disc radius was smallest. This may not be that 
surprising since streamlines of matter in the disc will be 
squeezed more closely together towards a region of smaller 
disc radius, causing extra viscous dissipation. We may ten- 
tatively conclude that cooling proceeds on a timescale less 
than the orbital period of the matter in the disc, which is a 
few minutes. 



4 DISCUSSION 

4.1 A central spike in AM CVn? 

The mysterious central spike has so far been observed in 
long-period, low mass transfer rate AM CVn stars, where 
the accreting white dwarf dominates the optical fiux. It is 
thought to originate on (or very close to) the surface of the 
accreting white dwarf, since the spike perfectly tracks the 
expected movement of the accreting white dwarf relative to 
the bright spot in the accretion disc. In addition, there is 
no accretion disc component that is expected to move at 
such low velocities, and the phase of the spike relative to 
the bright spot does not match the phase of the secondary. 
See iMorales-Rueda et al.l ll2003l) for a detailed analysis of 
the central spike in GP Com. 

The kinematic feature observed in the He I 4471 line in 
AM CVn, as presented here, has all the characteristics for 
being the perfect analogy t o the central spike in GP Com 
jMorales-Rueda et al.ll200^ . V396 Hya (Ste eghs et al. in 
preparation) and SDSS J124058.03-015919.2 jRoelofs et alJ 
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Figure 9. The accretion disc as a function of superhump phase. Left: the measured bright spot velocities in 16 superhump phase bins. 
Middle: the measured Hell 4686 emission centroids, again in 16 phase bins. Right: the bright spot velocity amplitude, together with the 
angle (in spatial coordinates) between the bright spot and the He II 4686 emission centroid, as a function of superhump phase. In this 
last panel we assume approximately stream velocities in the bright spot and Keplerian velocities for the He II 4686 emission. 
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Figure 10. Trailed spectra and linear back-projection Doppler tomograms of the Hell 4686 line as a function of superhump phase. One 
full superhump cycle is shown from left to right. 



12005^ : (1) it is intrinsically redshifted; (2) it agrees in ampli- 
tude with the expected velocity amplitude of the accreting 
white dwarf (and not with any other binary component's 
velocity amplitude); and (3) it agrees in phase with the ex- 
pected phase of the accreting white dwarf, relative to the 
bright spot. We therefore conclude that it is very hkely that 
it is the central spike. 

The mass ratio for AM CVn implied by the central spike 
is significantly higher than previou sl y tho ught, although 
there has been a study by IPearsonI i2003l) suggesting, in 
fact, a mass ratio g '--^ 0.19 — 0.25 based on a model in which 
the secondary star is moderately magnetic. The motivation 
for that study was to try to reconcile the mass ratio with the 
A'- velocity of the He ll 4686 emission line measu red (but dis- 
carded) bv lNelemans. Steeghs fc Grood (|2Q^), which sug- 
gested a mass ratio in this regime. It is clear from this 
work however, in particular figure EH that the measured K- 
velocities of the Hell 4686 emission indeed cannot be used 
reliably as an indicator for the motion of the accretor, since 
the apparent kinematic signal could easily be dominated by 
the eccentricity of the disc rather than the actual motion of 
the primary. 



4.2 Implications for AM CVn's formation channel 

One of the long-standing questions regarding the AM CVn 
stars is how they are formed. Theoretically, they can be 
formed in both double-degenerate and single-degenerate 
configurations: Roche-lobe overflow from a white dwarf onto 
another (more massive) white dwarf (e.g. iNelemans et alJ 
l ioflli), or from a heliu m-burning star onto a white dwarf 
Jiben fc Tutukovlll99J) . A third suggested formation chan- 
nel is Roche-lobe overflo w fr om an e volved mai n-sequence 
star onto a white dwarf ijPod siadlows ki. Han fc Rappaporl 
|2003l), where mass transfer starts by the time hydrogen core 
burning ends. This latter scenario, however, has the prob- 
lem of leaving significant amounts of hydrogen, while our 
high-S/N average and phase-binned spectra show exclusively 
helium and metals. For the other two formation channels, 
the main discriminator is the mass of the secondary star: 
in the helium-star channel, one expects a relatively hot and 
massive, only partially degenerate secondary, while in the 
double-degenerate channel, one expects a 'cold' degenerate 
donor star. 

The significantly higher mass ratio found for AM CVn 
in this paper has important implications for the nature of 
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the donor star, and hence for AM CVn's formation chan- 
nel, as it suggests a relatively massive donor. If we com- 
bine the orbital period with the mass-radius relation for a 
relatively 'cold' (core temperat ure Tc ^ 10^ K)^ degen erate 
helium white dwa rf as used in iNelemans et all ll200lll and 
refined recently in lPelove. Bildsten fc NelemansI (l200.'^^ . we 
get a secondary mass M2 ~ 0.035Mq and hence, through 
our minimum mass ratio of q = 0.18±0.01, a maximum 
primary mass of only M\ = 0.19Mq. This is very low; the 
maximum combined mass of M\ ~\- M2 ~ 0.22Mq would, 
for instance, be significantly lower than the average value 
Ml + M2 = O.SMr;^ found in detached WD-WD binaries 
fe.g. INelemans et al.ll20()5ll . This casts doubt on the possi- 
ble white dwarf (WD) nature of the donor star, and thereby 
on the double-degenerate formation channel for AM CVn, 
unless the white dwarf donor was still very hot upon Roche- 
lobe overflow - that is, unless mass transfer started shortly 
after the second common-envelope phase. 

Concurrent with the measurement of AM CVn's mass 
ratio presented here, it was discovered that AM CVn's dis- 
tance is larger than expected, through an HST/FGS par- 
allax measurement (tt = 1.65 ± 0.30 mas, Roelofs et al. in 
preparation). This means a smaller absolute magnitude My 
and thus, since we expect the accretion disc to dominate the 
optical flux, a higher accretion luminosity. 

In order to link the observed My to the mass of 
the secondary M2, we proc e ed in the same way as in 
iDelove. Bildsten fc Nelemanj i2005l) . We assume conserva- 
tive mass transfer that is driven entirely by loss of angular 
momentum due to the emission of gravitational waves and, 
since accretion proceeds via an accretion disc, we assume 
that all angular momentum carried by the transferred mat- 
ter is fed back to the orbit. We then have 



Parameter 



Value 



Ma _ J 2 

M2 ^ J C2 + 5/3 - 2q 



(5) 



where J is the orbital angular momentum, J is the orbital 
angular momentum loss rate 



J 
J 
the 



32 MiM2(Mi +M2) 



(6) 



due to the emissi on of gravitational waves 
l|Landau fc Lifsch'it^ll971^ ■ and 



C2 = 



d log R2 



(7) 



d log M2 

We ta ke C2 = -0.06 for M2 < O.2M0 as in INelemans et alJ 
i200ll) , based on evolutiona ry calculations for a mass- 
losing helium-burning star bv iTutukov fc Fedoroval il989f) . 
Clearly C,2 = Q around the point where helium burn- 
ing stops, while it may be somewhat closer to ("2 = 
—0.19 along the semi-degenerate part of the track for 
a secondary that is more helium deplet ed, as found by 
ISavoniie. de Kool fc van den HeuveJ lll986ll . This represents 
a <10% uncertainty in eq. (|5J|. 

The last step is to calculate the bolometric luminosity 
L, for which we use 



L = _M2($(Li) 



$(J?i)) 



(8) 



where $ is the common Roche potential, at the inner La- 
grange point Li and the surface of the accretor This 
states that the matter in the disc stays virialised during the 



Porb (s) 1028.7322 ± 0.0003 ('Skmman et al. 19991 

q 0.18 ±0.01 

Mi(Mq) 0.68 ±0.06 

M2(Mq) 0.125 ±0.012 

M2(MQ/yr) 6.71^;^ • lO-^ 

i (degrees) 43 ± 2 



d (pc) 



606'' 



(Roelofs et al. in prep.) 



Table 4. The collection of system parameters for AM CVn. 



accretion process, which should hold for approximately Ke- 
plerian particle orbits. 

Figure [TTj shows My as a function of M2 for q = 0.18, 
C,2 ~ —0.06 and a bolometric correction BC = —3.0 (see 
below). The lower limit to M2 is given by the require- 
ment that the system does not eclipse. The data point 
Mv = 5.0 ± 0.4 as determined from our aforementioned 
HST parallax is overplotted, which yields a secondary mass 
M2 = 0.125 ± 0.012 Mq and, via q = 0.18 ± 0.01, a primary 
mass Ml — 0.68 ± 0.06 M©. The corresponding mass trans- 
fer rate is M2 = 6.7li | • lO~^M0/yr, and the inclination of 
the binary comes out at i = 43±2°. Table|l]lists the system 
parameters for AM CVn in convenient tabular form. 

Our results have been corrected for the geometric effect 
that the 'apparent absolute mag nitude' ^IWarne^l^95^ of an 
accretion disc depends on the inclination. A factor 



-2.5 log 



cos I 
T/2 



(9) 



has been taken into account as the correction from apparent 
absolute magnitude to absolute magnitude, where the factor 
1/2 represents the direct ion- averaged fraction of the disc 
area seen by the observer, (cos i). We neglect the (unknown) 
effect of limb darkening. 

In converting the bolometric luminosity Afboi ob- 
tained from equation © to an absolute visual magni- 
tude Mv, the bolometric correction is of some impor- 
tance. AM CVn's spectral energy distribution from the far- 
UV to the optical ( Roelofs et al. in preparation; see also 
iNasser. Solheim fc S cmionofI 2001) indicates that it is dom- 
inated by a blackbody of ~30, 000 K. This corresponds to 
a BC of —3.0, with an estimated error of 0.3. If we calcu- 
late the minimum temperature the accretion disc must have 
in order to be able to radiate away the required luminosity 
(for any given Mi, 2 we know the maximum radiating sur- 
face of the disc), we find Tdisc ^ 30, 000 K for accretion rates 
M2 ^ 4 • lO"^M0/yr. The BC we use is thus self-consistent 
with the accretion rate we derive; in particular, it seems un- 
likely that we are overestimating the mass transfer rate due 
to an overestimate of the BC. 

Summarising, we have two independent pieces of evi- 
dence suggesting a relatively massive donor star - the large 
mass ratio and the large luminosity. If we compare our 
results to evolutionary models for helium star and white 
dwarf secondaries as given by.Nclcmans et aL ( 200 1») , we see 
that our derived secondary mass is compatible with the sec- 
ondary being the semi-degenerate core of a formerly helium 
burning star, while a 'cold' degenerate donor star is ruled 
out. See figure WIl The secondary mass M2 — 0.125Mq is 
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Figure 11. Relation between M2 and My for AM CVn. Our 
measurement of tlie absolute magnitude My and its implications 
for the secondary mass M2 are indicated by the dashed line. The 
open circle represents the minimum secondary mass for which 
the system does not eclipse (assuming a disc that extends to the 
3 : 1 resonance radius); the dotted line, towards larger M2 and 
smaller My, shows the helium-burning regime where we expect 
our assumptions regarding the mass-radius relation to break down 
(>10% error). The labelled dots indicate logM2 in Mq/jv. 



close to the turn-off m ass for helium burning, which lies at 
M2 ~ 0.16M^ iSavoniie de Kool fc van den Hc uvcl 1 98^ 
iTutukov fc FedorovafM89f) . We conclude that it is very likely 
that AM CVn formed through an evolutionary channel in 
which the donor star was, upon the start of Roche-lobe over- 
flow, still burning helium. These observations thus provide 
the first observational evidence that the helium star evolu- 
tionary channel contributes to the AM CVn population. 



4.3 Implications for AM CVn's gravitational wave 
signal 

We expect the orbit of AM CVn to be circular, so 
that the gravitational wave polarisation amplitudes A+ 
and Ax at twice the orbital fre quency, /, become (e.g. 
iTimpano. Rubbo fc Cornishll200d) 



Ay = -4 



c*d 



1 + cos i] 



(10) 
(11) 



where M = {MiM2f^^ /{Mi +M2)^/^ is the so-called chirp 
mass, i is the inclination of the binary, and d its distance. 
Defining the strain amplitude 



h^iliAl + Al 



1/2 



(12) 



log(M2/MQ) 

Figure 12. Relation between M2 and R2 for AM CVn. The up- 
per and lower solid line represent the evolutionary tracks for semi- 
degenerate helium star and for cold white dwarf secondaries, re- 
spectively <Nelemans et al. 2001: Tutukov & Fedorova 1989). The 
dotted line is the requirement that the secondary star fills its 
Roche lobe. The dashed line represents our measurement of M2. 
The data favour a semi-degenerate helium star secondary, rather 
than a white dwarf. 



and converting to more convenient solar units, one finds 



h = 2.84 ■ 10"^^ Vcos* i + 6 cos2 i + 1 



M 

Mq 



5/3 



^orb 

Thr 



-2/3 



1 kpc 



(13) 



Filling in the numbers derived above gives h - 



9 1 +0.4 
^•^-0.3 ■ 



10" 



Figure lT^ shows the estimated strain amplitudes and 
frequencies of the 'known' AM CVn stars. AM CVn itself 
is the first potential LISA source which is sufficiently well- 
constrained that we can plot error bars with some confi- 
dence, and we see that it stands out significantly above 
both the instrumental design sensitivity and the estimated 
confusion-limited Galactic background due to, mainly, de- 
tached WD- WD binaries. 



4.4 The superhump excess— mass ratio relation 

It is understood from numerical simulations, and it is ob- 
served in outbursting dwarf novae and nova-likes, that there 
exists a relation between the superhump period excess e and 
the mass ratio q of the binary. See Patterson et al. (200^, 
and references therein, for an overview including the very 
latest empirical work on the e(g) relation. In short, they 
conclude that for hydrogen-rich systems. 



e{q) = 0.18(7 -f0.29g^ 



(14) 



For AM C Vn, this results in g = 0.10 based on the pho- 
tometry of ISkillman et all l)l999|) . Although it is unknown 
wha t the intrinsic spread i n q{t) is, it appears from figure 
^2005^ to be quite low, no more than 



wha t the intrinsic sp : 
9 in lPatterson et al] 
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bg f (Hz) 

Figure 13. Gravitational-wave strain amplitude of AM CVn 
(solid circle). The upper and lower dashed lines show the de- 
sign sensitivities of LISA for a S/N of 5 and 1, respectively, in 
one year of data-collecting iLarson c t al. 2000). The solid line is a 
population syn thesis prediction for t he confusion-limited Galac- 
tic background iNelemans et al.l2004h . The open circles represent 
'known' AM CVn stars that are as yet not well-constrained in 
terms of cither distance or component masses, or both. 



about 0.01. This means that the g(e) for AM CVn devi- 
ates significantly from our minimum kinematic mass ratio 
g = 0.18 ±0.01. 

Of course, the superhump excess e may depend on more 
parameters than just the mass ratio q; it is quite conceivable 
that the helium nature of AM CVn may be of influence. A 
recent study by iGoo dchild fc Qgilvid 1120061) indicates, for 
instance, quite a significant dependency of e on the thick- 
ness of the accretion disc, where thicker discs lead to smaller 
superhump period excesses for a given mass ratio. If this ef- 
fect should be at play, our results would suggest a physically 
thicker disc in AM CVn compared to discs in hydrogen-rich 
systems of the same mass ratio. 

Although AM CVn represents only a single data point, 
it seems advisable to apply the empirical superhump excess- 
mass ratio relation with caution for such ultra-compact bi- 
naries. An obvious test would be to obtain the mass ratio 
in HP Librae, the other known helium nova-like, and check 
whether it, too, shows a discrepancy with relation 11411 . 



4.5 Size and eccentricity of the accretion disc 

In our derivation of the minimum mass ratio q — 0.18±0.01, 
we have assumed a maximum accretion disc radius equal to 
the minimum radius of the primary Roche lobe. It is custom- 
ary to assume a slightly more stringent limit on the maxi- 
mum accreti on disc radius, caused by tidal truncation of the 
disc. See e.g. I Warn er (1995) and references therein. For ref- 
erence, we show in figure ITU the often-used tidal truncation 
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Figure 14. Same as the left panel in figure |S] but now including 
the tidal truncation radius (eq. 1151 '> for an inviscid disc, indicated 
by the dot-dashed line. 
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Rt 
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0.6 



(0.03 < g < 1.0) 



(15) 



for an inviscid flow, which is indeed smaller than the full 
minimum Roche lobe radius we assume. However, since this 
tidal truncation radius depends, in particular, on the un- 
known viscosity of the matter in the disc (and increasing 
with viscosity) , we have not used it as a hard upper limit to 
the disc's radius. Interestingly, one sees that the best-fltting 
accretion disc radius we derive for ballistic stream velicities 
in the bright spot (reproduced in figure lies very close 
to, but slightly above, the tidal truncation radius HSU . 

The effective eccentricity value e = 0.04±0.01 we derive 
for the outermost part of the disc, from the stream-disc 
impact spot, is smaller than typical values e ~ 0.1 — 0.2 
found in numerical simulations. Although the eccentricity is 
not one of the observables that is usually quoted in numerical 
studies (the resulting superhump period excess is often the 
only given quantity), it can be inferred that the effective 
ec centricity we find is signifi cantly smaller than that found 
bv lSimpson fc WoodI il998l) in their simulations of discs in 
AM CVn stars. 

A possible explanation would be that tidal truncation, 
discussed in the previous paragraph, acts to circularise the 
disc at the outer edge to some extent by effectively re- 
moving particles from the most eccentric orbits there. This 
tidal truncation process is essentially a competition be- 
tween viscosity transporting angular momentum outwards 
and tidal forces 'dissipating' this angular momentum in- 
creasingly strongly with radius, and predictions about this 
tidal truncation process would thus be affected by our lack 
of understanding of the viscosity in accretion discs. 



12 G.H.A. Roelofs et al. 



4.6 Spin of the accretor and tidal synchronisation 



5 ACKNOWLEDGMENTS 



If the central spike originates from the surface of the accret- 
ing white dwarf, in particular if it rotates with the accretor, 
its width puts an interesting constraint upon the accretor's 
spin. The central spike is unresolved in our data, and we 
therefore place an upper limit of ~45 km/s, equal to our 
formal spectral resolution, on the FWHM of the feature. 
For the system parameters derived for AM CVn, corotation 
with the orbit would imply a projected equatorial velocity 
of ~36 km/s. The central spike we observe is thus consistent 
with a corotating primary, while significantly faster rotation 
is ruled out. 

We can estimate the tidal synchronisation timescale of 
the accretor from this constraint. For simplicity, we assume 
that AM CVn has been steadily accreting matter at a rate 
M2 ~ 7 ■ 10~^ AIq I yr for a long time. As long as the accretor 
is near corotation with the orbit, the spin-up timescale r 
of the primary due to the accretion of disc matter can be 
written as 

T = — f5i — k\ — : (16) 

to Povh \ G M2 

where k ^ 0.18 fsee iMarsh. Nelemans fc Steeehl |2004) is 
the moment of inertia factor [k = 2/5 for a solid sphere), 
and u! is the angular velocity of the accretor. Requiring that 
tidal synchronisation proceeds on timescales similar to or 
shorter than the accretion-induced spin-up of the primary, 
gives a maximum tidal synchronisation timescale Ts < 2 • 10^ 

yr- 

A key question in binary evolution theory is whether 
two white dwarfs, upon Roche-lobe overflow at orbital pe- 
riods of a few minutes, can manage to stabilise the ac- 
cretion process and avoid a merger by feeding back an- 
gular momentum from the spun-up accretor to the orbit. 
This has profound implications for the number of systems 
that may survive the initial phase of mass transfer, when 
accretion proceeds at a high rate of M2 ~ 10~^MQ/yr 
jMarsh. Nelemans fc Steeghsll2004 . 

Although its absolute magnitude is highly uncertain 
from a theoretical point of view, tidal synchronisation is un- 
derstood to scale with the orbital separation and the radius 
of the accretor as 



(17) 



,M2j \Ri, 

(see, e.g., iMarsh. Nelemans fc Steeehsl l20o3 ) . This means 
that, for a three- minute white dwarf binary in which Roche- 
lobe overflow is about to commence, the tidal synchronisa- 
tion timescale will be at least a factor 10'' shorter than for 
AM CVn, or Ts < 200 yr. This is an interesting result since 
it is exactly the regime that is needed for an appreciable 
fraction of binary white dwarf systems to survive the ini- 
tial ph ase of mass transfer. See[Marsh, Nelemans & Stc eghsl 
i2004ri . in particular their flgure 11. 

If the central spike is thus associated with the ac- 
creting white dwarf as we think, it implies a tidal syn- 
chronisation timescale that is short enough to make the 
doub le-degenerate formatio n channel for AM CVn stars vi- 
able jNelemans et aljl200j) . It would thereby also increase 
the viability of the direct-impact accretor scenario for the 
ultra-compact binaries RX J0806. 3-1-1527 and V407 Vul 
jMarsh fc Steeghsll2002h . 
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